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Abstract 

We present a novel and rigorous approach to the Langevin dynamics of ideal 
polymer chains subject to internal distance constraints. The permanent con- 
straints are modelled by harmonic potentials in the limit when the strength 
of the potential approaches infinity (hard crosslinks). The crosslinks are as- 
sumed to exist between arbitrary pairs of monomers. Formally exact expres- 
sions for the resolvent and spectral density matrix of the system are derived. 
To illustrate the method we study the diffusional behavior of monomers in 
the vicinity of a single crosslink within the framework of the Rouse model. 
The same problem has been studied previously by Warner (J. Phys. C: Solid 
State Phys. 14, 4985, (1981)) on the basis of Lagrangian multipliers. Here 

we derive the full, hence exact, solution to the problem. 
61.41. +e, 64.60.Cn, 87.15.By 
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I. INTRODUCTION 



A theoretical treatment of the dynamics of polymer networks is a generally unsolved 
problem. In a preliminary attempt Edwards et al. [jl]J^ studied the problem of a polymer 
subject to internal distance constraints. In their investigation the underlying theoretical 
problem was to handle the quenched degrees of freedom (hard crosslink constraints) which, 
for example, in a random network exist between pairs of arbitrary polymer segments (the 
monomers). As a first step Edwards considered a (macroscopically) long polymer chain 
which was internally crosslinked to itself at random. The polymer backbone was assumed to 
be Gaussian and the resulting dynamics was found to be of the standard Rouseian type 
Permanent junction points were treated by Lagrangian multipliers, which led to enormous 
technical difficulties for the corresponding differential equations. In fact, these could only 
be handled by strong approximations, such as pre-averaging in combination with harmonic 
variations. Even when the problem was highly oversimplified and only one crosslink was 
considered the method of Lagrangian multipliers still becomes highly involved as was pointed 
out in a successive paper by Warner 0. 

The purpose of the present paper is to develop an alternative formalism for treating 
Langevin dynamics of polymers subject to internal distance constraints. For calculational 
simplicity the simplest working model for a free polymer, the Rouse model I^Q, is considered. 
It is suggested that the more complicated problem of a random network can also be treated 
exactly by the presented method. We adopt here the minimal model suggested by Edwards 
imj^ and consider one (macroscopically) huge polymer chain which is randomly crosslinked 
to itself. Such a crosslinking process will lead to tetrafunctional crosslinks. In previous works 
we have already demonstrated that the analogous static problem can be solved exactly when 
excluded volume effects between the polymer segments are ignored Physical quantities 
such as the static structure factor or the radius of gyration were found to be self averaging and 
could be determined by relatively simple numerical means. The essential trick was to account 
for the crosslinks in a general connectivity matrix, that includes both the connectedness of 
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the polymer chain and an additional contribution from the crosslinking. In close analogy 
we expect the corresponding dynamic problem to have a similar exact solution as long as 
complicating factors such as excluded volume, hydrodynamic forces or entanglements are 
neglected. To demonstrate this analogy we start from the standard Langevin description for 
the polymer segments and solve the stochastic differential equation in terms of its resolvent. 
As an instructive example we reconsider the eight-shaped polymer problem (i.e., a polymer 
ring with one crosslink) studied by Warner and present its full solution. We first confirm 
the results from Warner, which have been derived only for low frequencies and low Rouse 
mode index, but show secondly the exact solution in the entire frequency and mode domain. 
Moreover the technique introduced here opens new ways to study Langevin dynamics of 
constrained systems. 

The paper is organized as follows. In section II the physical model - a generalized version 
of the Rouse model with internal distance constraints - is introduced. Section III summarizes 
some of the basic theorems regarding Langevin dynamics to be used later on. In section IV - 
the main calculational body of the paper - the general mathematical formalism for handling 
internal distance constraints is developed in detail. Our treatment is a generalization of a 
method previously developed for computing statistical properties of randomly crosslinked 
Gaussian structures, i.e., ideal polymer networks In section V an application of the 

method to diffusional motion of a single crosslink is given (the Warner problem). Section 
VI contains a short discussion of main results and outlook. 



II. ROUSE MODEL WITH INTERNAL DISTANCE CONSTRAINTS 

As a minimal model for the dynamics of a Gaussian chain subject to internal distance 
constraints we consider a generalized version of the classical Rouse model [0. Its discrete 
version is a bead-spring model, where the motion of the beads (monomers) is governed by 
the coupled set of Langevin equations 
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In this equation of motion the inertial term is omitted as usuaL ( denotes the inverse 
mobihty or friction constant, and Ri(t) {i = 0, ...,N) are the trajectories of the monomers 
in 3-dimensional space. The stochastic forces Fj(t) are assumed to be 5-correlated with first 
and second moments given by |^ 

(Fnt)) = , (2) 
(t')) = 2CKTS,,S^f3Sit - t') . 

Superscripts a, P = x,y,z represent the 3-dimensional Cartesian coordinates. In the classical 
Rouse model excluded volume interaction and hydrodynamic forces are disregarded and only 
elastic forces between monomers are retained in the Hamiltonian. Here we consider a more 
general form of the Rouse model with an extra potential to allow for modelling the internal 
distance constraints 

o N o M 

/?^o = E(R^ - R-i)' + - • (3) 

1=1 e=l 

The first term in the Hamiltonian represents the connectivity of a Gaussian chain with 
persistence length a, whereas the second term models the crosslinks. In particular we are 
concerned with permanent constraints when a monomer, say ii, is linked to another monomer 
labeled by ji . For more than one crosslink a whole set C of crosslink " coordinates" is needed 
to specify all junctions in the system 

C = jl), (ieje), {iMju) ■ (4) 

For example, depending on C the object under investigation can be a flexible ring polymer, 
a two-dimensional membrane, or a rubber network (figure 1). 

Although the theory will be developed for arbitrary coupling constant e two scenarios 
are of special relevance. For e — > it has been shown that the Hamiltonian (3) is suitable 
to model hard (5-constraints (the classical crosslinks) of the form 

M 

nw.(t)-R>(t)). (5) 

e=l 
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The case e ^ oo leads to the well-known problem of a free chain which serves here as a 
reference state. One might be worried that the above model is ill-defined and might diverge 
in the limit e ^ 0. It will be shown in section IV that the converse is true and that a 
surprisingly simple solution can be obtained for this special limit. As shown in the earlier 
paper on the static properties, ref. |]^, it is important to take the limit e ^ at the very 
end of the calculaton. This procedure ensures firstly, that no mathematical problems occure 
and secondly that in this case hard crosslink constraints is treated properly. Before going 
into more of the calculational details some of the basic definitions and notations regarding 
Langevin dynamics are summarized in the next section. 



III. LANGEVIN DYNAMICS OF IDEAL POLYMERS: PRELIMINARIES 

Consider the generalized Ornstein-Uhlenbeck process specified by equations (l)-(3). 
For calculational simplicity matrix notation will be used. We define + 1 dimensional 
"super- vectors" with three dimensional vector components to account for the positions of 
all monomers R(t) = (Ro(t), RAr(t))^ and for the stochastic forces acting upon them 
F(t) = (Fo(t),...,F7v(t))"^. The dagger denotes the conjugate complex of the transposed 
vector. Furthermore, the + 1 dimensional connectivity (Kirchhoff) matrix is introduced 
as 



1 M 
^ e=l 



(6) 



where 



Wo 



1-1 ■■■ 

-12-1 ; 

■•. ■•. ■•. 

; -12-1 

■■■ 0-1 1 



(7) 



is the Wiener matrix associated with the polymer "backbone" , and 
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(8) 



models a single crosslink. For further use we note that a characteristic time scale is given 
by the inverse of the "frequency" 

3KT 



(9) 



The dimensionless parameter z = {s/aY in equation (6) is used to enforce the crosslinking 
constraints. With the above definitions the system of stochastic differential equations (1) 
cast into matrix form reads 



^ + A^(.)R(t) = lF(t). 



(10) 



Some of the physical quantities of interest and their interrelations are listed below. More 
details can be found for example in reference 0. The Green's function to equation (10) is 
given by 



g{t) = lime--^(")* 



(11) 



For z — > the case of hard (5-constraints is recovered. Otherwise z is an additional distance 
parameter in the model. Of great importance in the following derivation is the Laplace 
transform (resolvent) of the matrix A4 



roo , \ — 1 

7^(cu) = / dt e"*'^* Q{t) = lim (iul + M {z) 

Jo z^O V 



(12) 



where I denotes the identity matrix. 



/^From equation (2) the spectral matrix of the Langevin forces F(t) is found to be 

{F{uj)F^{uj')) = Uirk^TC 6{uj ~uj')I , (13) 

where 

/oo 
dte-'^'F{t) (14) 
"OO 

is the Fourier transform of the stochastic forces. By FjFj we mean the usual 3-dimensional 
scalar vector product, whereas FF^ is used for outer vector products. Fourier transforms 
are always denoted by tilde. 

A formal solution to equation (10) can be obtained by Rice's method f^. The spectral 
density matrix for the stochastic variable R(t) can be derived by use of (13) and Fourier 
transformation of (10) 

(R(w)R^(cu')) = 127tD S{u - uj') n{uj)n\uj') , (15) 

with the diffusion coefficient D given by 

D = k,T/C . (16) 

Of primary interest for the diffusional behavior is the two-time correlation function matrix 
defined as 

C{t,t') = ((R(t) - R(0)(R(t) - W)^) ■ (17) 

Finally a steady-state solution for C{t, t') in terms of the resolvent (12) is easily derived from 

the expression for the spectral density matrix in (15) 

ion r<xi . , 
C{t - i!) = / dLo(l- coscu(t - t')) n{uj)TZ{uj)\ (18) 

TT JO ^ ' 

In the following study our primary goal will be to find a general approach to calculate the 
resolvent 7?.(u;), equation (12), for an arbitrary set of crosslinking constraints C, equation 
(4). From there Green's function and correlation functions can in principle be obtained by 
use of the standard formulas presented in this section. Although M.{z) is a matrix which 
highly depends on all the details of C (the crosslink positions), substantial progress can be 
made by invoking the following exact method. 
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IV. CALCULATION OF THE RESOLVENT n{u) 



The first step in deriving a general expression for TZ{u!) for hard crosshnks is to find 
a way to perform the hmit ^ — > in equation (12). This is an interesting problem in its 
own right which so far could only been handled by introducing a finite cutoff at z = 1 and 
successive crude variational estimates. Here we present an analytically exact approach that 
can overcome these difficulties. The mathematical trick is to utilize an additional symmetry 
of the crosslink term in (8) by writing the complete crosslink contribution in (6) in form of 
a dyadic (outer vector) product 

M 

^A'(ie,Je)=WW^ , (19) 

e=l 

where 

W(C) = (ui,...,Um) (20) 

has been introduced as the {N + 1) x M rectangular matrix with each of its M column 
vectors given by 



, (e=l,...,M). (21) 



Here represents the + 1 dimensional unit vector with 1 in the ieth position, and 
otherwise. Thus U{C) has only 2M elements not equal to zero that contain complete 
information about all crosslink positions. In the above notation each crosslink is uniquely 
represented by a vector Ug. Note that all vectors Ug, e = 1, M are linearly independent 
for tetrafunctional crosslinks. Combining equations (6), (12) and (19) the resolvent cast in 
matrix form reads 

n(uj)^—lim(i—I + yVo + -UU^)~\ (22) 

It is convenient to decompose (22) into a singular and nonsingular part with the nonsingular 
part being 
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W = i—I + Wo. (23) 

It is well-known in the mathematical literature that if the inverse of W exists, then the 
inverse in equation (22) is given by 



(W + ^UU^) ^ = W-^ (24) 



z 

X 



This theorem can be directly verified by matrix multiplication. The latter identity is also 
known as Sherman-Morrison formula M. 



A. The limit z ^0 

There are two subtle points about the existence of the right hand side of (24). First we 
require W^^ to exist. The only critical case arises if u; = 0, i.e., when W = Wo in (23). The 
problem here is that Wq is only positive semi definite and there is one mode with eigenvalue 
from translational invariance. This can be directly seen from the definition of Wo in (7) 
which is a row (column) constant matrix. However, even in the semi definite case the above 
theorem remains valid if W~^ denotes a generalized inverse of W as was proved in reference 

Secondly from the definition of Ug in (21) it is easily verified that for tetrafunctional 
crosslinks all M vectors Ug are linearly independent. Thus in general the kernel zX+U'^W^^U 
will be a positive definite matrix of dimension M and full rank which has only positive 
eigenvalues for all nonnegative values of z. As a consequence performing the z — > limit in 
equation (24) leads to a well-defined expression for the resolvent 

n{uj) = —w-Ui -uiu^w-^uy^u^w-^) . (25) 

The first term is the linear chain (Rouse) model, whereas the second part arises entirely 
from the effect of crosslinking. Although the case of general crosslinking potential z is still 
implicit in the basic formula (24), we will restrict ourselves in the following discussion to 

9 



the somewhat simpler case z = 0, i.e., hard (5-constraints. Equations (24) and (25) are 
formally exact solutions to the problem posed in equations (l)-(4). The further evaluation 
of TZ{u!) for specific realizations of crosslinks C can be split into two parts and is discussed 
in subsequent sections. 



B. Resolvent of Wq 



Evaluation of the inverse of W in (25) can in principle be done by full diagonalization 
of Wo which is tridiagonal. For calculational simplicity we consider here only the cyclic 
counterpart of Wo with periodic boundary conditions 



Wo 



2-1 ■■■ -1 
-1 2 -1 ■■■ 



V 





-1 



-1 2 -1 
0-1 2 



(26) 



Both models (7) and (26) are known to obey the same Rouse dynamics in the limit N oo 
||1 1|| . Physically the latter situation represents a flexible ring polymer. The eigensystem to 
(26) is of particular simple form since it is a circulant. The eigenvalues read 

irk 



Ai, = 4 sin 



N+1 



k = 0,...,N . 



(27) 



The modal matrix of (26) is the Fourier matrix JF ||l4l with matrix elements 



2nikl 



kl 



exp 



k,l = 0,...,N . 



(28) 



y/N+l + 1 

Spectral decomposition leads to the well-known representation of the inverse W~^ in terms 
of its eigenvalues 

1 



kl 



N pY-n 27rm(fc-0 



(29) 
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C. Discussion of kernel 



The remaining calculational task for determining Tl{uj) is the evaluation of the kernel 
function in the second part of (25) 

]C{lo; C) = {U^W-^Uy\ (30) 

Since /C(u;; C) depends on all the crosslink positions C=(zi, ji), {iM.iM) via U no further 
analytical progress is possible without specifying the crosslink in the system. On the other 
hand, from the mathematical structure of lC{uj] C) most problems of interest fall into one of 
the following three categories. Only one of these will be considered in detail in section V. 

(i) The number of crosslinks M is small. Since /C(u;; C) requires inversion of an M x M 
matrix analytical progress is always possible if M is not too large. A particularly simple 
problem is treated in the next section when we consider the dynamics of a polymer shaped 
like the figure-of-eight (figure la). 

(ii) Another special case arises when M is large, but there is some additional pattern in 
the structure of lA. Examples of this kind are illustrated in figures lb and Ic. In particular 
the sketch in lb shows an example of a macromolecule with distance constraints z 7^ 0, 
i.e., the more general case governed by equation (24). For the above examples /C can be 
calculated as a consequence of the regularity of the crosslink positions. We will report on 
these systems in a separate publication. 

(iii) The third important category arises when M is large and the crosslink positions C 
in (4) are picked at random. This is the case of a polymer gel (figure Id). Here one has to 
resort to numerical computation of K,{uj\ C) [0. However, there is still a huge calculational 
advantage with (25). For a polymer network we have in general M <^ N. Equation (25) re- 
quires "only" the inverse of an M x M matrix 0] and not of the complete N xN connectivity 
matrix as is commonly believed in the polymer literature [|12|,|13|. An analytic approach to 
the network problem would be to perform the quenched average of the resolvent over the 
crosslink positions C. The latter problem is a key problem in current network research and 
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has not been analytically solved even for the static problem. 

Before calculating TZ{uj) for a specific example, we want to establish some remarkable 
and general properties of the operators in equation (25). Consider the crosslink part in (25) 

V = w^^uiu^w-^uy^u^w-^ . (31) 

By elementary matrix multiplication it is found that 

VWV = V , U^V = U^W~^ , VU = W^U . (32) 

A matrix with these properties is said to be a generalized projector to W~^. Furthermore 

{WVf = WV , (VW)2 = VW (33) 

are idempotents whose eigenvalues are known exactly: Ai = 1 and A2 = with degeneracies 
M and N — M. By use of the above results it is easy to prove that the resolvent satisfies a 
remarkable orthogonality relation 

U^n{uj) =n{uj)U = . (34) 

Equation (34) is valid for arbitrary crosslink positions C and independent of the specific 
crosslink topology of the system under investigation. 

V. DIFFUSIONAL MOTION OF A SINGLE CROSSLINK 

As the simplest possible application of the method developed in section IV we consider 
the figure-of-eight shaped polymer depicted in figure la. What we have in mind is to model 
the dynamics of a single crosslink in an ideal dilute network when the distance between 
crosslinks is large . It is expected that monomers in the neighborhood of the crosslink are 
somewhat affected by the slower dynamics of the crosslink p|flp!5[] . 

A suitable realization of the system in figure la would be 
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U{C) = Ui = eo — ejv+i . (35) 

That is, monomer is linked to monomer (A^ + l)/2. The main calculational task is to 
determine the kernel function, equation (30), of the system. Prom (35) and with given 

by (29) we get immediately 

where the summation includes only the odd terms. For the diffusional motion the quantity of 
interest is the self-correlation function contained in the diagonal elements of the correlation 
matrix (18) 

[C{t-t')U = {{Rs{t)-Rs{t') f) , (37) 

where s is the distance of the sth monomer with respect to the crosslink at position s = 0. 
Typical terms and manipulations in the straight-forward derivation which is not carried out 
in detail are of the form 

1 , exp(2nisx) 
'ax- 



-I 

2 h 



iuj/ujQ-\- 4 sin {ttx) 



In deriving the first integral we have performed the N ^ oo limit. The latter expression 
was obtained by setting sin(7rx) ~ ttx. Only the final result for the self-correlation function 
(37) is quoted here 



[Rsit) - Rs{t')f) = A{s, \t - t'\) 2a'^^^^^-—- . (39) 



t-t' 



The time-dependent prefactor is given by 



A{s,t-t') = l =/ dx 40 



X e *'^(^ cos (s'Vx) sin (s'^x)) 

= l-.lrrfyOOs(jg!-^)'°S(\+y'). (41) 

An Jo V 2 y3/2 ' v ; 
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which scales with 



s' = 



2s2 

(42) 



The comphcated integral in (41) is plotted in figure 2. The asymptotic behavior for small 
values of s' is governed by the expansion 

1 S\/? S'2 



In particular for the diffusional motion of the crosslink (s = 0), we find A = 1/2 which 
is exactly half the diffusion constant of an unconstrained monomer in the Rouse model 
||16|| . The finding is in agreement with the result in reference based on the method of 
Lagrangian multipliers. For s —>■ oo (monomers that are sufficiently far from the crosslink) 
we recover the diffusion law of the classical Rouse model {A = 1) which was first derived by 



de Gennes 16 



In addition we obtain the crossover from the slower dynamics of the crosslink to that 
of a "free" monomer in the classical Rouse model as s is varied from zero to infinity. The 
crossover takes place on time scales of the order 

Ts = s^ujo = {safUm , (44) 

where sa measures the distance of the monomer from the crosslink (figure 3). Interestingly 
a monomer begins to feel the presence of the crosslink only after a timespan of the order r^. 

The two limiting cases A = 1/2 relevant for the slower dynamics of the crosslink and 
A = 1 for the "free" chain segments far away from the crosslink are expected on physical 
grounds (dashed lines in figure 3). An "inner" chain segment has only two neighbors, 
whereas the crosslink is surrounded by four neighbors. Thus in general for a monomer with 
functionality / a prefactor A{f) = f/2 is expected as was pointed out previously 

VI. CONCLUSION 



Within the framework of the Rouse model we have proofed that an exact solution for 
the Langevin dynamics of a polymer subject to hard delta constraints exists when excluded 
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volume and hydro dynamic forces are neglected. The fundamental and general result for the 
resolvent, equation (25), was derived for an arbitrary crosslink configuration making our 
results also applicable to the challenging problem of a random network. 

In this investigation we restricted ourselves to the simplest physical scenario, where only 
Rouseian dynamics was involved. This case was deliberately chosen to highlight the principal 
mathematical difficulties. As a special application we studied the dynamics of the figure-of- 
eight shaped polymer depicted in figure la. In contrast to an earlier attempt by Warner ||5| 
based on Lagrangian multipliers which yielded only two limiting cases s — > and oo, the 
full solution could be derived by our method. Moreover, our result allows for computation 



of the dynamic scattering function [T^ and comparison with experimental data taken in the 
dry network state. A detailled comparison will be studied in a future, less formal paper. 

For the physically more realistic scenario of a swollen network in a Theta solvent further 
generalizations are required like taking hydrodynamic interaction into account. A general- 
ized version of the equation of motion (1) would read 



= E (-Vr,.Ho({R,}) + F,(t)) , (45) 

where @ij is the Oseen tensor. Although the above equation becomes analytically un- 
tractable, for most experimental situations a pre-averaged treatment is well justified 0. 
This computation is left for future work. 
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FIGURES 

FIG. 1. Examples of different crosslinking topologies, (a) The polymer shape discussed in 
section V. a is the persistence length of the polymer backbone. The hard crosslink constraint is 
enforced by £ — >■ 0. (b) A ladder shaped polymer with e ^ 0. (c) Two-dimensional membrane, (d) 
Random network. 

FIG. 2. Prefactor A{s'), equation (41). The dashed lines represent the asymptotic behavior for 
small and large values of s'. 

FIG. 3. Crossover of the mean squared displacement from the dynamics of a "free" Rouseian 
monomer (upper dashed line) to the slower dynamics of a crosslink (lower dashed line) . 
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